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RADIATION OF CYLINDRICAL DUCT ACOUSTIC 
MODES WITH FLOW MISMATCH 


by S. D. Savkar 
with I. H. Edelfelt 

General Electric Corporate R&D 
SUMMARY 


Calculations for the radiation of spinning acoustic modes, 
with or without a centerbody, and with or without flow tempera- 
ture and velocity discontinuity, are presented. Solutions to 
the appropriate convected wave equations devised around Fourier 
transforms and Wiener-Hopf technique are presented. The decom- 
position of the asymmetric kernel, resulting from a flow and 
temperature mismatch, is carried out in part exactly and 
partially using the so-called Carrier-Koiter approximation 
procedure. The resulting solutions offer a good approximation 
to the radiation of both symmetric and asymmetric modes through 
a flow discontinuity represented as a plug flow jet issuing 
from a cylindrical duct. Besides the Koiter approximation, 
the major limitation on the calculation program is the diffi- 
culty of calculating the high order Bessel functions with 
sufficient accuracy. 

Sample calculations for symmetric modes are presented 
offering comparisons between the current calculations and in 
particular available experimental data such as that of Plumblee 
et al [6] and Mechel et al [7] . So long as one restricts 
attention to the amplitude functions, such as the radiation 
pattern or the absolute value of the conversion (reflection) 
coefficients, there is generally a qualitative agreement, even 
here there are certain peculiar features of the predictions 
which do not appear to be reproduced experimentally for reasons 
discussed. The agreement between the calculated and measured 
phase values in the presence of flow appear to be rather poor. 

In the case of the asymmetric modes, limited comparisons 
between calculated and measured radiation patterns appear to 
be reasonable, thus in general reinforcing the conclusion 
regarding ability to predict amplitude functions. However 
none of the data available go to significantly high enough Mach 
number. For this one must rely purely on calculations. 

The calculations suggest that the influence of the Mach 
number discontinuity is to offer a conflict between the effects 
of convection and refraction. The former becoming increasingly 
more significant at high subsonic Mach numbers by the presence 
of a strong lobe in the shadow zone. Effects of flight condi- 
tions are to weaken such conflicts so that in the limit of 
uniform Mach number, one recovers Carrierte well-known solution. 

The effect of a temperature mismatch is rather more 
dramatic, resulting in a sharp beaming of sound off to the 
side. This is consistent with the kinematic picture of a 
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plane wave impinging 


a temperature discontinuity. 


I . INTRODUCTION 


In their now classic paper, Tyler and Sofrin [1] demon- 
strated that turbomachinery noise propagates as spinning 
acoustic modes and is eventually radiated to the far field. 
While the qualitative features of the physical process were 
brought out, their use of a radiation field uncoupled from the 
process within the duct and no flow left a number of unanswered 
questions. There appeared to be a need for the solution of the 
complete radiation problem coupling the field within the duct 
to the far field in the manner of the Levine and Schwinger [2] 
paper and its extension to a duct immersed in a uniform flow 
by Carrier [3] . From the point of view of turbomachinery 
application, the need was to be able to handle high lobe num- 
ber spinning (asymmetric) duct wave guide modes and the flow 
discontinuity present as for example in the exhaust mode. The 
first poses computational difficulties since it requires an 
ability to handle high order Bessel functions accurately. 
Lansing et al . [4] have developed a theory for the asymmetric 

modes in extension of the Carrier paper [3] . However the use 
of an isothermal, uniformly flowing medium excludes such fea- 
tures as refraction and an enhanced mismatch between the duct 
and the surrounding. The refractive influence of the flow 
mismatch was demonstrated by Mani [5] , who employed a form of 
the Carrier-Koiter approximation scheme. Unfortunately these 
calculations are restricted to flat ducts. 

While the theory has been slow in evolving, the experi- 
mental measurements of relevant model problems have not been 
easy to come by. Since the problem of primary interest must 
include a flowing medium, the measurements have been limited 
by the inadequacy of acoustic instrumentation in turbulent 
flow. Then there is the problem of generating suitably well 
defined acoustic modes, especially the asymmetric spinning 
modes of interest. Perhaps the most comprehensive attempt 
along these lines is that of Plumblee et al . [6], who have 

developed an annular flow duct apparatus capable of handling 
flow velocities up to 200 ft /sec. together with the acoustic 
instrumentation necessary to map the acoustic fields both 
within the duct and in the radiation field. Plumblee et al's. 
[6] attempt at measuring spinning mode radiation impedances 
are unique. One of the very few other measurements available 
are those of Mechel et al. [7], who measured the plane mode 
reflection coefficients both with and without flow. When it 
comes to more realistic measurements, as from fans, the only 
available data are for the radiation field. Even these in 
general do not offer a good basis for discerning the influence 
of specific modes and moreover the picture is clouded over by 
the incoherent scattering of the noise by flow turbulence and 
velocity gradients. The problem is clearly more severe in 
the exhaust mode than the inlet mode. Crigler and Copeland 
[8] have made some rather detailed measurements of a fan 
radiation field which appears to be strongly dominated by a 
rotor-IGV interaction of a definitive modal character. 
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In this report we shall attempt a solution to the problem 
of radiation of spinning duct acoustic modes of high lobe order 
and include in our consideration a flow mismatch both in the 
sense of a jet exhaust flow and a temperature mismatch between 
the jet fluid and the surrounding. The problem in a sense re- 
presents the fan exhaust or turbine exhaust problem. However 
we shall model this using a cylinder of zero wall thickness 
from which issues a plug flow jet through which upstream sound 
refracts and radiates to the far field. We shall also approxi- 
mate the influence of a centerbody by allowing the inner 
boundary to be a doubly infinite "rod". The solution will be 
devised around the appropriate convected wave equation using 
the particle displacement matching boundary condition and the 
Wiener-Hopf technique. In the process of the solution, we 
shall resort to the Carrier-Koiter approximation scheme to 
render the problem of decomposing the Wiener-Hopf kernel 
tractable. The validity of the scheme has been extensively 
discussed and illustrated by both Koiter [9] and Carrier [10] . 

One notes also the successful use of the same by Mani [5] in 
the flat duct problem. 

We shall attempt an evaluation of our calculation routine 
by comparing sample predictions with the earlier discussed and 
other experimental results. We expect at the outset reasonable 
correlation for the radiation patterns but not for the reflection 
coefficients, or more precisely the phase of the reflection 
coefficient due to difficulties resulting from experimental 
conditions and more importantly the presence of turbulent flow. 

Finally, we shall also offer some interesting results 
related to temperature mismatch. 


NOMENCLATURE 


a 


H^ 1 ^ ( Z ) 
m 

I (Z) 
m 


J (Z) 
m v ' 


K m< Z > 

M 

P 


Y ( Z ) 
m 


radius of duct, used to nondimensionalize all 
the length parameters 

speed of sound 

Hankel Function of the first kind 
modified Bessel Function of the first kind 
/^1 

Bessel Function of first kind order m 
reduced frequency, oja/c 

modified Bessel Function of the second kind 

Mach number 

pressure 

Bessel Function of the second kind 
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a Fourier transform variable 

a hub to tip ratio of the centerbody 

p density 

to circular frequency 

<J> velocity potential 

y mg s x zero of the appropriate Bessel function or cross 

product order m 

n particle displacement 

II. THEORY FOR THE RADIATION OF SPINNING DUCT WAVE GUIDE MODES 
II. 1 The Exhaust Mode - Formulation 


Consider now the problem of a monochromatic spinning 
acoustic mode lobe number 'm' in semi-infinite duct of radius 
'a = 1.' The duct wall is assumed to be of a negligible thick- 
ness and is located at r = 1 . , x < 0 as in Figure 1, where 
(r, 6, x) represent a cylindrical coordinate frame. The ana- 
lysis will be restricted to small disturbances in perfect, 
inviscid, nonheat conducting media. In the formal analysis, 
however, an allowance will be made for a mismatch between the 
flow within the duct and the surrounding medium. The region 
formed by the duct and its phantom extension will be designated 
as region 'I' or '1' in the subscripts. The surrounding region 
is then II or 2 in the subscripts. In what follows the fre- 
quency of the incident signal, propagating in the duct toward 
the open end, will be prescribed to be 'w' and the time de- 
pendence of the form exp(-joot), j = /^l, will be assumed. 

Then the disturbance velocity potential amplitude cf> is governed 
by the convected form of the Helmholtz equation 

V 2 <f> + (k + jM ^) 2 4> = 0 (1) 

where it is understood that M = in region I and M = M 2 in 
region II and that k = k^ = uia/c^ and k = k„ = u)a/c 2 respectively. 
The latter makes an allowance for differences in temperatures 
(or gases) between regions I and II. The solution to the above 
problem is to be carried out subject to the boundary conditions: 

a) Along the walls of the duct the transverse acoustic 
particle velocity be zero, i.e. 

=0 r = 1 and x < 0 

b) The acoustic pressure be continuous on the extension 
of the duct into x > 0, 
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M 1 3(f> l M 2 3<f> 2 

+ j kj lx ^2 + J “a# ’ r = 1 - 

c) In addition to the above, noting the definition 

Dr] _ 

Dt ar * 

it is required that the radial acoustic particle 
displacements in the two regions be matched along 
the streamlines separating the two regions. This is 
the boundary condition proposed by Ribner [12] and 
Miles [13] . There are some interesting questions 
related to this boundary condition, but we defer 
these to a later section. 


Adopt now a representation in terms of the Fourier trans- 
forms, so that in the two regions I, II, we write, taking a 
circumferential dependence of the form exp(j m 0): 


and 


*1 


j A + (o) 


T m (Y l r) 
Y 1 


e -J«x 


da 



oo 

B + (a) 

J 

— CO 


V^2 r) 

*2 K m<V 


e -jax 


da 


( 2 ) 


(3) 


where a is the Fourier transform parameter I m and being 
modified Bessel functions of the 1st and 2nd kind 

and Y 2 =a 2 -(k+aM) 2 (4) 

It is easily verified that equations 2 and 4 and 3 and 4 
satisfy the wave equation. The '+' signs in (2) and (3) 
designate functions regular and free of zeros in the upper 
half a plane, to be defined presently. In the procedure we 
have adopted, we shall draw heavily on the book by Noble [11] . 
Hence to simplify repeated references, such as to specific 
equations or pages, we shall adopt the notation (N, 1.15) in 
reference to equation (1.15) of Noble or page (N-10) as 
referring to page 10 of Noble [11] . 

Now the function y of a clearly has two branch points at 
k/(l - M) and -k/(l + M), so we introduce branch cuts such 
that y -*• +°° as a -> ±°°. As discussed by Noble, the cuts can 
be made in several different ways, however for our subsequent 
calculations the cuts which turn out to be convenient are 
illustrated in Figure 2, where for the sake of illustration 
it is assumed that k has a small, but positive imaginary 
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component, i.e.* 

k = k r + i e , taking a = £ + j r) . 

The cuts are along the hyperbolas defined by £n = k e , in the 
limit e -*■ 0 these reduce to cuts such as illustrateS, for 
example, by Figure 3. Actually, it turns out that the cuts 
for y-^ could just as conveniently have been made to extend 
horizontally from k^/Cl - M^) to +°° and -k^/Cl + ) to 

but we adopt the above convention for consistency. This allows 
us now to define the upper and lower half planes (taking M ^ , 

M 2 to be positive) as 

upper half plane - for Im(a) > the larger of Im(-k 2 /(l + 

M 2 )) or Im(-k 1 /(l + M 1 )) 


and 


lower half plane - for Im(a) < the lesser of Im(k 2 /(1 - 

M 2 )) or Im(k 1 / ( 1 - M^) 

Consequently the path of integration in (1) and (2) is taken 
to pass between the branch cuts in the region of commonality 
(shaded region in Figure 2). 

Clearly the form of solution adopted ensures the satis- 
faction of the boundary condition (a). Hence consider an 
incident wave of the form 


Pine = V^mn r > ex P<J K mn x > 


(5a) 


or 


<j> . 

r mc 




J (y r ) 
m mn 7 


“ o - M 1 >W k 1 > 


exp(j K x) 


mn 


(5b) 


where K mn the axial propagation constant is defined as 


♦Writing a = £ + j n , the branch is defined as 

9 1 + 0 2 
J( ••• o - 


with the ranges 


Y = I Y I e 

-2tt + arc tan A + < 0^ < arc tan A + 


arc tan A < 0g < arc tan A + 2ir 


where 

.+ = n - £/(l - M) 
C - k r /(l - M) 


- = n + c / ( 1 + M) 
K + k r /(l + M) 


A 


7 


K 


ran 



a - MMn 

' — / r~ 


mn 




( 6 ) 


with the branch chosen being Im(K m ) •> 0, corresponding to 
downstream propagation, y beingthe zero of J'(y ) = 0, 
satisfying boundary condition (a). Then the ra9iaT n acoustic 
particle displacement is 


or 


n • 

me 


n = 


j 9 4>/ 9 r 


o)( 1 - M K /k) 
v mn' 7 


mn 
w^(l - 


J ’ ( y 
m v mn 


r) 


M. K 


mn 




exp(j K x) 


mn 


(7) 


By its definition the incident wave particle displacement 
vanishes as r + 1_ for all x. Again, the acoustic particle 
displacements corresponding to (2) and (3) are 


= 1 

0) 


| A + (a' 


I ' (Yi 
m ' 1 


r) 


(1 + a M 1 /k 1 )I r ;(Y 1 ) 


-jax , 
e J da 


( 8 ) 


= 1 

0) 


B + (a) 


k ;^ 2 


r) 


(1 + a M 2 /k 2 )K;( Y2 ) 


-jax , 
e J da 


( 9 ) 


Substituting into the boundary condition (c) and noting that 
hinc E 0 for r = 1_, there obtains the relation between A + and 

B + 

A + (0t) B + (a) 

(1 + a J^/k^ ~ (1 + a M 2 /k 2 ) = F - (a) (10) 

The left and the right hand sides being regular and zero free 
respectively in the upper and lower half planes. Noting the 
behavior of the two sides as x -> 0, by the standard application 
of the Wiener-Hopf argument and the Liouville theorem, each 
side is equal to a constant taken to be zero. There obtains 
the relation 


(1 + a M /k 1 ) 

A + (a) " (a + a M 2 /k 2 ) B + (a) 


(ID 


The remaining boundary condition, continuity of acoustic 
pressure 


p 2 (i + > - + P lnc (iJ 


for x > 0 


noting the definitions 
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oo 

r 


Pi = 


J P 1 W 


(1 + a M- /k. )I (y. r') 

V a) v Wv r 1 1 e JOX da (12) 


Y 1 *ii<V 


(1 + a M /k )K (r r) 

P 2 J p 2 “ J B + (a) Y 2 K^(y 2 ) 6 dC * ^ 13> 


and on r = 1 for x > 0, writing 


me 


■> J m (p mn ) 
2 tt 


f e~ Jax 

J < a + K mn> 


da 


(14) 


since 


j K x 
mn 


= _ _X_ f .jC 

2tt j J a + 


-jax 


K 


da 


mn 


then yields the functional equation 


p 9 (1 + a M 2 /k 2 )K m (y 2 ) (1 + a MJk,) 2 I m ( Yl ) 


l'~l 


m v ' 1' 


^2 K m<V 


(1 + a M 2 /k 2 ) Yl r( Yl ) 


x B + (a) 


J (y ) 
m v mn 


271 Pl a.(a + K mn ) 


= G (a 


(15) 


where the equality is understood to apply to the region of 
commonality (shaded region of Figure 2), and the function 
G_(a) is regular and free of zeros in the lower half a plane. 
Denote by L(a) the grouping in the curly brackets. Then 
equation (15) may be rewritten as 


L(a) B + (a) 


J (v ) 

m v mn 


2tt p 1 w(a + K mn ) 


= G (a) 


(16) 


Then if the function L(a) can be decomposed into factors 
L + (a) and L_(a) regular and free of zeros in the upper and 
lower half planes 


L(a ) = L + (a) L_(a) (17) 

The standard application of the Wiener-Hopf arguments and the 
Liouville theorem yields 


R _ m v mn 

+ v J 2t r p. to(a + K )L^(a) L (-K ) 

1 mn' + v 7 - mn' 

Except for the decomposition of L(a), this completes the 
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formal solution to the problem posed. 


II . 2 Decomposition of the Wiener-Hopf Kernel 
II. 2. a The Carrier-Koiter Approximation 

Recall the definition of the kernel function L(a): 

T , , p 2 (1 + “ VWV <1 + « W 2 W 

' P 1 y 2 W ' Cl + » M 2 /k 2 ) Yl I'( Yl ) 


Formally speaking the required split can be carried as per 
theorem c, page (N-15), of Noble, whence 


*n(L + (a)) = -gfj 

— oo 

with the integration along the 
hatched region in Figure 2 


&n(L(Q) 

?-o 


lower boundary of the cross- 


and £n(L (a)) = r 

the integration path being along the upper boundary of the 
crosshatched region of Figure 2. However, it does not appear 
possible to carry out the integrations explicitly, or reduce 
them to simpler integrals as would be the case if L(a) were 
purely an even function of a. In view of this, the decom- 
position was carried out using the Carrier-Koiter approximation 
scheme [9, 10]. The applicability of this approximation to 
the acoustic radiation from a flat duct has been previously 
demonstrated by Mani [5]. To motivate the approximation, 
regroup L(a) as follows. 


+ jr 


+ jr 


ftn(L(g)) 
5 - a 


dC 


L(a) = 


where 


(1 + a M 1 /k 1 ) 2 2 

(i + a M 2 /k 2 ) Yl i;(y 1 )y 2 k;( Y2 


) 


- V 


T(a) 


(19a) 


( y — y ) 

T(a) = 2e 2 1 y 2 I m (Y 1 ) K^(y 2 ) H(a) 


(19b) 


H(a) 



(1 + a M 2 /k 2 ) 2 
(1 + a M 1 /k 1 ) z 



X 


w w A 

W " J 


(19c) 
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We will presently need the limit Mm H(a), which is readily 
shown to be a °° 


Mm H(a) 

a oo 



(19d) 


assuming the same static pressure within and without the jet, 

P 1 c3 = p c* and where and M are the ratio's of specific 
heats in regions II and Irrespectively. In the case of uniform 
flow, = M 2 and P 2 /P 1 = 1, the square bracketed quantity in 
(19a) is essentially a normalized Wronskian. Consider the 
substitute kernel L*(ot) defined as 


L*(o) = 


(1 + o M 1 /k 1 ) 


(1 ♦ C M 2 /k 2 ) y x i;(Y 1 )Y 2 k;(y 2 ) 
with T*(a) being defined as the quotient approximation 


P " 1 

4 


" V 


) 


T*(a) (20a) 


T*(ot) 


6 B + a + 6/g 
X a 7 + a + 6/ X 


(20b) 


With — = Mm [-/yg/Y^ H(a)] it is clear that L*(a) + L(a) 

a -*■ 00 

in the limit as a ->• 00 along the real axis. A more complete 
match according to the Carrier-Koiter rules is then obtained 
by matching L(a) and L*(a) at a = 0 up to its first moment: 


6 = T( 0 ) 


* n _ T(eo) - T( 0 ) 
6/A T ' ( 0 ) 


(21a) 

(21b) 


«/3 " ( 6 M)f££) (21c) 

The objective now is to carry out the decomposition of L*(a). 
To facilitate the subsequent discussion, group L*(a) as 


L*(a) 


1 M°0 f 4 (« 

8 ? *(,) f 2 (a) 


where T*(a) is defined by equation (20b) and the functions M 
are defined as 

Y-, -m + 1 -y 

f l (a) = ( ^" ) e 

Y 9 m + 1 y 

f 2 (c0 = e K m (Y 2 ) 

f 3 (a) = (1 + a M 1 /k 1 ) J /(l + a M 2 / k 2 > 
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Y 9 in 

f 4 (a) = (^) 

In terms of this simplified notation, the integral decomposition 
of L*(a) is readily carried out. For an example, L*(a) is 
simply expressed as 


£n[L*( 




00 + js 


- 00 + js 


An[f 3 (c)] 
l - a 


dC 


00 + js 


-OO + js 


An[f 4 (c)] 

? - a 


dC 


• r dS *n f l(C ) 


? - a 


dC 


-oo + js 


* f jS *n[f 2 ( C) ] 


-OO + js 


C - a 


dC 


One can similarly write down the factor L*(a). Each integral 
will now be treated separately in the following subsections. 


II. 2. b Decomposition of T*(a) and fg(a) 


The decomposition of the quadratic quotient T*(a) is 
trivial, the roots of the quadratic terms being 


R 1) 2 ~ ~ \ ~ ' /;L / 4 ~ S/B (22a) 

R 3 , 4 = - | ± /1/4 - 6/A (22b) 


Quotient groupings such as for example 
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a - R a - R 

/ ^ 7 X f^tc or / ^ T etc - 

3 4 

being assigned to the appropriate half planes by testing the 
location of the roots R. . 

l 

The decomposition of the remainder kernel can be carried 
out exactly. Since the factor (1 + a M 1 /k 1 ) 2 /(l + a M 2 /k 2 ) 
or f 3 (a) is clearly a plus function, hence write 

fg(a) = f 3 (a) 

f“(a) = 1. 

II. 2. c Integral Decomposition of f 2 (a) 


The grouping f 2 (oi) was motivated by the suggestion made 
by Noble (page N-130). It should be noticed that f 2 (a) is in 
fact a function exclusively of Y2> hence we may write 


f 2 (°0 


f(Y 2 ) = 



Y < 


K m<V 


(23) 


The function f(Y 2 ) is free of zeros in the entire complex a 
plane. Integral decomposition of this function can be further 
simplified by making the following transformations in the a 
plane 


a 02 1 M 2^ 

k 2 = k 2 /(l - M 2 ) 4 

in terms of the Sq 2 plane, f(y 2 ) an even function of 5 q 2 . 
This at once allows us to simplify the integrations involved 
by using the procedure 'c' of Noble, page (N-19). Hence deform 
the contour of integration in the a 02 plane as per Figure 3, 
then 

+ g (5 Q2 ) 

f 2 = f + = {f(Y 2 )> 4 e 2 (25) 

-g< a Q2 ) 

f 2 = f_ = (f(Y 2 )} 4 e 2 (26) 


k„ M, 


(1 - M 2 ) 1 


(24) 
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V»e omit here the details of the algebra involved in deducing 
an expression for g(0 with simply a reference to Noble [11] 
(pages N-19 and N-20). Then for |£| £ £ 2 , it is easily verified 
that 


g(0 = 


25 p [ 

* 6 << £ 2 ' 


fl(u)u du 


K 2 ) - u 2 } /k 2 - u 2 


P 


u du 


0 


{(k 2 


? 2 ) - u 2 }/k 2 - u 2 


ii f ^3 + y2) 
tt I y 2 + K 2 
o 


dy 


A 2 + y 2 

tt j y 2 + P 

o 


dy 


with 


- ~ ^ Z) + i¥r ^2 ■ 

fi(x) = arc tan(Y^(x)/J^(x) ) - ^ 


(27) 

(28) 


'P' preceeding the integrals denoting that they are to be taken 
in the sense of principal value. The integrations can be con- 
siderably simplified and carried out explicitly for k 2 suf- 
ficiently greater than the first zero y ml of J^(x) = 0, by 
using the asymptotic expansions of ft(x) given by Abramovitz and 
Stegun [14] , thus 

/ \ ,m , l x , 4m 2 + 3 , 16m 4 + 184m 2 - 63 + ... 

Sl(x) * x - ( 5 + 3) + - 8x ■ * 384? 


Usefulness of the above expansions is well illustrated by the 
sample comparisons for m = 0 and m = 10 in Figures 4 and 5. 

It appears to suffice that k 2 be more than roughly halfway 
between y m ^ and y m2 . Hence obtain for 0 £ E, < k 2 and k 2 suf- 
ficiently > y m i, using only the first 3 terms of the above 
expansion 


g(5) = 

it k 0 


0 


sin0 ^(k^ sin0) - /l - £ 2 /k 2 B(/k 2 - £ 2 ) 
(sin 2 0 + £ 2 /k£ - 1) 


d0 


14 


- _ . 2 4m 2 + 3 

s «3 _ 


r TT 

to - COS 


17 2 


- 1 #-i 


(29) 


Provided k 2 is sufficiently larger than y m i , one may readily 
deduce other similar formulae. We shall present here only two 
additional formulae of interest in deducing the conversion 
coefficients in mismatched flow, thus first for £ > k„ > y r 
we have 


ml ‘ 


k. 


g(0 


= _ 21 

TT 


ft(x)x - k 2 + 


{(x 2 - k| + £ 2 )/k 2 - x 2 


dx + £ 


- ¥ [(1 + V - % + ?, ** (#- * /Fjq 

17 2 4 2 8/T z - k 2 k 0 2 


- 1)] (30) 


and for £ 2 = k 2 > y 2 

2 H m, 1 


g(5) = f 


fl(x)x - X 2 


0 


{(k| - C 2 ) - x 2 }v^| - 


dx 


+ jr- (| + + f Mi_L_3] 

2 4 8 k 2 


(31) 


II. 2. d Product Decomposition of f^(oc) 
Consider next the factor f^(a) 


-Y- 


f 1 (a) = 




<Y,/2) 


m - 1 


-Y- 


(32) 


The factor e can in fact be extracted and treated separately. 
Using the branch cuts illustrated schematically in Figures 2, 3, 
the decomposition of e can be carried out in a manner par- 
allel to the procedure used for f 2 (ot). Since the results are 
identical, for the sake of brevity, we quote instead results 
adopted from Noble (N-1.35) making the transformation 


a 


01 


= a ( 1 - M 2 ) 4 - 


k l M 1 

(1 - M 2 )^ 


k ± = k 1 /(l - M 2 )* 


(33) 
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Thus the exponential function regular and free of zeros in the 
upper half plane for 0 i “01 - k l is 


exp{ — arc cos(a n1 /k 1 )} 


01' “1 


(34a) 


with its analytic continuation for 

Y 1 S oi + Y 1 
exp{- j -± *n(-^ i)} 

11 k 

K 1 

so that in the limit ® in the upper half plane, 

behaves as 


(34b) 


it 


a 


01 


exp{-j -jp ln(2 d ()1 /k 1 )} 


(34c) 


Jhere now remains the factor ^ r J 1 ('Y^)/(Y^/2) 
f^(ct) its inverse 


m - 1 


^(a) 


(Yi/2) 


m - 1 




(Yi/2) 

JyTT 


m - 1 


m 


Denote by 


(35) 


recalling that y. = -j y 1 for -k./CL + M^) < a < k., /(I - MA 
In this form the^decomposit ion or J-,(a) is readily accomplished 
using the infinite product form of J^CY-i)- We shall again 
omit the details of this well documented procedure with re- 
ference to Noble, pages (N-128) and (N-129), and simply note 
the results: 


~+ |2(Y 1 /2) m N - a 

ft(a) =\ ± ir 

' 1 J i<V s = 1 a i, s + a 


where 


a 


01 


+ j Y 1 (a ) 


(36a) 


. Ski - y 7 (1 - Mf) - k- M, 

_ 1 H m, s 1 ' 1 1 

“l, s (1 - MJ) 


a. 


- K. s < l - M i> + k i M i 


■1, s “ (1 - M£) 

N being the highest number for which 6L are real, and 

X , s 
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V“> - - a -ir 11 - c + *" r + i - F^-)' 


s — 1 


m, s 


00 a 

1 


1 Y 

+ l arc sin 


a 


01 


s = 1 M m, s s = 1 


Z 2 


<°i. s + 5 W> 


with 


a* = u 2 - k 2 
1 , s m , s 1 


and 


C = Euler constant = 0.5772 ... 


One can write down a similar expression for f (a) or simply 
note that 1 

fl(a) = ij(a) fj(a) (36b) 

The function Y^(a) can be put into a form more convenient for 
computation by expanding the arc sin into a series and using 
the Digamma function, ip(Z). Thus, T(Z) here being the Gamma 
function , 

<KZ) = r ' (z)/r(z) 


and 


00 

1 A- s4-z - c + *<i ♦ z> 


s = 1 


find 
Y 1 (a) = 


a 


^ [1 + iKN + | + |) + £n 2ir/k 1 ] 


a ni 00 

01 I 


IT 


it u , . XT , m 3 

s = Is + N + t- t 

2 4 


✓a 


^ + Qt Z 

1, s + N a 01 


^ 1 A 2r + 1 a 01 { 

r = 1 s=l/a z 


2r + 1 
■) ) 


+ a; 


1, s + N ^01 


with 


A 2r + i = t 77 (2s - 1 ) ] / [ 2 r r!(2r + 1)] 


s = 1 


( 37 ) 


The choice of the particular form of the grouping now becomes 
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clear. Since in the limit a + °° in the upper half plane, for 
example, the factor regular and free of zeros in the upper 
half tends to (Noble page (N-128)). 


r- -Y- 


fj(a) = 




(Yj/2) 


m - 1 


^01 ^ 

exp[-{x(a) T-p l n 


(1 - C) + (-^ + 5 - f + ^)£n(-^i + ? + t) 
J tt j tt 2 4 2 ' v jir 2 4 


a 


01 




J1T s = 1 S %, s 


Hence together with the remainder of the kernel, the trailing 
edge condition is satisfied by picking x (°0 to be 


a 


X(a) = 


01 


j a 


01 


[1 - C + In (P-) 


+ 


s 




(38) 


II. 2. e Decomposition of f^(a) 


This now leaves just one factor to be decomposed f (a) = 
( Y 2 / 2 ) m / ( Y^/2 ) m which was introducted in conjunction with the 
earlier factors. 


and note that 


*4(«> 

* 4 ( a ) 


« , >. / T 2 .m 

f 4 (°o = (— ) 

Y i 


= m( 


Y 2 


n. 

Y 1 ‘ 


Hence obtain the factor f^(a) closing the contour in the upper 
half plane 
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f + 4 (a) = 


(V 

]2 


+ *02 

*1 

• 

" 5 0l " 


1 

J 2 

a 02 

*1 

+ 5 01 _ 


(39) 


with f 4 (a ) = f 4 (a)/f 4 (a). 

This completes the discussion of the Wiener-Hopf decom- 
position of L*(a) into two factors L$(a) and L*(a). These 
now allow us to readily deduce quantities of physical interest. 


II . 3 The Radiation Pattern and Conversion Coefficients 
II. 3. a Asymptotic Solution to the Radiation Pattern 

The simplest procedure by which to deduce the radiation 
pattern is to asymptotically evaluate the integral pressure 
expression in the limit (r 2 + x 2 )4 -»■ °°. The expression for 
the acoustic pressure field in region II obtains from (13) and 
(18) 


r J (y )(1 + a M n /k« ) 

_ • m mn /v 2' 2 

p 2 3 P 2 I 2tt p i cj(a + K )L*(a) L*(-K ) 

1 v mn ; + - mn 


K m ( ^2 r) 

K ; ( ^> 


-jax J 
e da 


(40a) 


It is convenient now to use the spherical coordinate frame 
(R, © , $) centered at the duct termination, then with 
reference to Figure 6 


x = R cos © 
r = R sin © 

Thus in the limit R -*■ °° (or more precisely R sin @ ->• «>) one 
can use the asymptotic expression for R sin (H) ) 

_ ^ -Y« R sin © 

K m (Y 2 R sin©) - A/2y 2 R sin © e 

There obtains the far field pressure approximation 


j p 2 f 

2 " p l J (o + 


Wn )(1 + “ M o/ k o> 


K mn )L*(a)LM-K mn )K'(Y 2 )y- 


~ 2 j 2 v ^/2 R sin W * 
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-j R h(a) 
e ° v ' da 


(40b) 


where h(a) is the free field propagation factor 

h(a) = a cos (§) - j y 2 sin®. 


(40c) 


The inversion of (40) can now be carried out by the method of 
stationary phase [15] . The point of stationary phase cor- 
responds to the root of the equation 

h ' (a) = 0 

straightforward algebra yields the requisite root a 


0 


a 0 (1 - M z ) [M 2 


COS 


a. 




sin z (fi) 


-] 


There obtains then the solution to the radiation field 


J (u )D ( (5) ) 
m VM mn mn v ^ J 


-j R h(a Q ) 


D mn<® ) being the directivity factor defined by 


mn 


j P 2 (1 + “o W^o + K mn>~‘ ^ 3/2 e-J ”/ 4 

2 p l L ^- K ™m>l h 0 ‘Vl* L i<VW 


(41) 


(42a) 


(42b) 


and where after some algebra 


h " (a ) = - 
0 ( 0 J 


(1 - M 2 sin 2 ® ) 3/2 


k 2 sin z (g) 


(43) 


The — tt / 4 factor arising in the directivity factor from the 

-j k sin ® 


negative sign on h^' (a^) 


and 


Y 0 = 


/I - si 


'I 


in" © 


The kernel functions L*(a) and L*(a) being 

V 


N 

TT 


*1, s - a 


* ^1 y 2 (y 2 ) s - i 


(44) 


+ a 
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fioi »2 *< 5 02 > . Y 1 . 

T +(“) ex P<-2 2 2 + V aro cos <“ 01 /k 1 > 


L*(a ) = 


r i(“) - 


+ j Y 

k 2 " a 02 k l + a 


01 


k 2 + & 02 k l " S 0U 


bi/2 


(45a) 




" Y i Y 2 J i<Y 1 )H 1 [ 1 1), (T 2 ) S = 1 a 


NS- + a 
it -U_§ x 


1 , s 


- a 


T*(a ) exp ( - 


“01 Y 2 . g ( S 02 ) . Y 1 


+ ^ + — (" “ arc cos(a 01 /k 1 >) 


- j V“> - 


(45b) 


Y-( ) and g(a Q2 ) are defined by equations (37), (27) to 
(31): Then supposing for example , R 3 are roots in the 
upper half a plane 


T*(a) = /W/^ 


a - R^ 
a - R, 


a - R, 


and T*(ct) = /Wf\ 


II. 3. b Evaluation of the Conversion Coefficients 

The conversion coefficients are deduced from the inversion 
of 


2ti 


(1 + a M- /k - ) 2 J (p )L* _1 (-K ) I (Y, r ) 

1 1 m VM mn y - mn' m w l e ~jax da (4g) 


(1 + a M 2 /k 2 ) 


(a + K mn )L»(a) 


y ! : ;< y i> 


with the closure contour in the upper half plane. There are 
no branch points, only the simple zeros of Yi I m Y ^i)- Denote 
these zeros, obtained by solving the equation Yj_ = j p , by 
to obtain for the pressure field within the duct in tne 
limit x -*■ 


p„ = J (p ) 
m mn y 


J m (>j m n r) J K mn x 
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s = 1 
th 


_ J (u r) -i K x 
n m in ins “ ms 

R — r— 7 r— e 

ns J ( u ) 
m ms 


where R" s is the m mode conversion coefficient of the n 
radial mode into the 


(47) 

th 


R 


m 

'ns 


(1 + 


K M 1 /k 1 ) ' 
ms 1' 1 


L* ‘(-K, 


)L| ‘(K„ 


v-1 


+ K ms V k 2> 


(K + 
v ms 


mn ' — + -ms )(1 " m ' v m S y 

) (k- M. - K (1 - M*)) 
mn v 1 1 ms v 1 


K 


(48a) 


and 


M, 


K 




ms 


(1 - 


M? ( 1 - 

ms 

IT! 


■mJT 


(48b) 


Kjjjg being the s x reflected mode propagation constant, with 
the plus sign on the radical being picked for the upper half 
plane . 


II . 4 Influence of a Doubly Infinite Centerbody 

Suppose now that the duct inner wall is a doubly infinite 
cylinder radius a, then adopt for p^ and p 2 the representations 


p l 


= J P-l lo 


A + ( a ) ( 1 + a M 2 /k 2 ) x 


V^L r > K r^l - K m^l a > -jax , 

V^V^l a)] 6 


(49) 


= j P, 


(j 


B 


+ (a) 


( 1 + a 


V k 2> K m^2 


r ) 


-jax 


k ;<v 


da (50) 


Before proceeding further, let us introduce some short- 
hand notation to simplify otherwise bulky expressions involving 
cross-products of Bessel functions which we shall repeatedly 
encounter in the annular duct calculations. Specifically, we 
standardize on the notation of equations (9.1.32) of reference 
[14]. However, since several of the symbols have been used 
previously in a different context, we shall use asterisks to 
indicate the difference. In addition, to avoid ambiguity we 
will always specifically indicate the arguments involved. 

Hence note the following definitions: 


<lm (a ’ 


b) = 


J (a)Y'(b) - J’(b)Y (a) 
m v ' nr ' m 7 m 


(51a) 
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r*(a, b) 
m v ’ 


s*(a, b) 


J’(a)Y (b) 
m v m v y 


J' (a)Y' (b) 
m m 


J m (b)Y'(a) 


J' (b)Y’ (a) 
m y m 


(51b) 

(51c) 


Using the above shorthand notation, the incident wave can be 
expressed as 


p. =q*(u r ,U a )e 
*inc M m mn ran 


j K° x 
J mn 


(52) 


Then repeating the earlier steps, the Wiener-Hopf kernel L(a) 
may be rewritten in the form 


L( a) = [L(a ) ] 




p 2 Y 1 (1 + “ V k 2> 




P 1 Y 2 WW (1 + ° W 


2 

T ~ 



^2 K ^V W 


(1 + a M 2 /k 2 ) 2 
(1 + a M 1 /k 1 )^ 


- 1 


(53) 


The first square bracketed term is precisely the kernel L(a) 
decomposed for the cylindrical duct to which reference may be 
made. There remains the last factor which we designate as 
W(a) for convenience. Again adopt the Koiter approximation 
W*( a ) 


W* ( a ) ~ 4r ~ ~ ~ - 4^~ (54 

A a 2 + a + 6 / A 

where by matching the function in the limits a 0 and a -*■ 00 
and the first derivative as a ->■ 0 we obtain 


and 


W(a ) = 1 

A a °° 

6 = W(0) 

6/X = w (°°) - W(0) 


( 55 ) 


W’ (0) 
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6/6 = 6/A 


W(0) 

Wl'oo'l 

"V / 


Now repeating the earlier calculations for the radiation 
field and conversion coefficients obtain the result, in the 
far field for p 2 (a Q being defined by (41)) 


P, 


q*(y 

^m mn 


mn 


a )D° (©) 


-jRh(a Q ) 

R 


(56a) 


with the directivity factor D° ((H)) 

mn 


D (©) = 

mn vxy/ 


£ p 2 (1 + a o V k 2» tt o + O ' 1 ^i 3/2 e ~ J,,/4 
2 Pl L ^ ( - K ran ) l h " (“o^ 4 ^<“o )K m <y 2 )/5TF ®~" 


(56b) 


where y now are the zeros of 
mn 

5 m s £ (0 ' 5 m " 0 < 56c > 

y - = 0 m f 0 are trivial and are not included in the set 
ox ’ admissible eigenvalues. 


To determine the pressure field within the duct, invert 
(49) with a contour closed in the upper half a plane. Again 
there are no branch points and the inversion is carried out by 
summing the residues at the zeros of Yi [ Im( Yi )Km(a Yi) - 
SiOjHiO Yi)]. Denote the values of a at these points by K a 
Then using tne cross-product notation of [14], as defined abov6, s 
the pressure field within the duct is given, away from the exit 
plane of the duct, by 


= q*( y 
m mn 


y o ) 

mn 


q* ( y 
M m v p mn 


mn 


o) j K, 


mn 


q*(y 
^m mn 


mn 


o) 


+ l 

s = 1 


.m 


R .i. , ^m^ms r ’ ^ms l ^ *Sns x 

ns \ q*(n , n a) }e 
^m ms ms 


(57) 


where the conversion coefficients R m are now defined as 

ns 


_ (1 + K° M 1 /k ) )L*“ 1 (K u )(1 - m 2 /y 2 ) -1 

r>m _ ms V l J - v mn' + v ms ' v /h ms' 

ft "" ■ X 

nS (1 + M 0 /k 0 ) (K° + K a )(k. M 1 - K a ( 1 - M 2 ) ) 

ms 2' 2 ms mn 1 1 ms v l y/ 


( y - m/a) r*(y " n a) 
i , _ ms ' ' morns’ M ms 2 

1 + o 7-=-z 


(y 2 - m z ) 

v ms 


q*(y , y a) 
rn v ms’ ms 


( 58 ) 
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II. 5 The inlet Problem 


Strictly speaking, the plug flow jet approximation is 
totally incorrect in the inlet mode. However allowing the 
license, one can treat such a formulation as a crude first 
approximation, especially in view of the intractability of 
the problem with a sink flow. For the inlet mode, the flow 
duct is considered to be located at r = 1 and x > 0. Again, 
both M 1 and M 2 are taken to be positive. 

The inlet problem now is very similar to the exhaust 
problem in that the governing equation and boundary conditions 
are quite similar to those already considered. Hence only an 
outline of the solution is given below. Recalling that 


_ I i£ = [JL + Me —1 2 n 
p 9r l 9t 9x J n 


Let 


n l = 


A_(a) I i ;(Yi r)e JaX dCX 


(59) 


(60) 


a M 


°° A (a)(l + 


Pi = 


P 1 W 


I )2 e -j“x 


V Y i r)da 


Y l I i <Y i > 


(61) 


(It is readily verified that (60) and (61) satisfy (59).) 


^inc 


exp(-j K* x)y J ' (y mri r) 
mn mn m mn 

{ M. K* ^ 

^ , 1 mnV2 

~^~J 


where 


K* = 
mn 


fk 2 . - (1 - M^)y 2 + k 1 M. 

1 v 1 mn 1 1 


(1 - Mj) 


(62) 


(63) 


and 


p. = exp(-j K* x)J ( y r) 
me ^ mn m mn 


(64) 


In region 1 1 : 


B (ct)K' ( y 0 
- v y m v ' 2 


r)e- JOX dc. 


K' (Y„) 
m v ' 2 


( 65 ) 
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Po = 


P 2 a) 


B ( a ) 1 1 + 


0 O 2 

* 1 e ” JUA K m ("Y 2 r ) da 


[ Mr, a"l 


*2 K m ( V 


( 66 ) 


As before (62) indicates that n- = 0 for r = 1 for all x. 
Hence by (65), (60), the Wiener-Hopf technique and Liouville's 
theorem : 


A_(a) = B_(a) (67) 

Application of the pressure boundary condition for x < 0 
yields : 


a M„ 2 


p 2 (1 + ~k 


0 z a M 2 

A K„(y,) (1 + -c-^) I m (Y x ) 


m v '2 


p l ^2 K A 2 > 


'll 


A (a) 




+ 2tt ( a - k£) Pi ^ - H + ( “> 
Again if we let 


L(a) = 


“ a M_ 2 

p 2 (l + -^) K m ( Y2 ) 

p l y 2 K m ( V 


(1 + 


a M„ 2 


kT> W 


'll W 


(68) 


(69) 


and if L(a) can be factored as L + (a) L_(a), then the solution 
of (68) is that 


A (a) = 


-j J (Vi ) 
m mn 


Pi 2 *<“ - K S„> E + ' K S„> E - ( “> 


(70) 


(70) completes the formal solution to the inlet problem. 

Note from (69) and (19) that 

a M 

L(a)/(1 + -7-4) = L(a) (71) 

2 

Since (1 + a M 9 /k 2 ) is clearly a plus function we have 
immediately that 
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L _(a) = L_ ( a ) 

and \ (72) 

L + (a) = (1 + a M 2 /k 2 )L + (a) J 

Equation (71) shows that any formula or method to_factorize 
L(a) will immediately yield the factorization of L(a). 

We wish to note that for the inlet problem if there is a 
uniform flow over the entire plane the application of the 
boundary condition that Ap(x) = 0 for x < 0, r = 1 (where Ap 
denotes the jump in acoustic pressure across the phantom duct 
extension described by r = 1 and x < 0) automatically ensures 
that Acf> is also zero across such a surface. This is seen by 
noting that if we write 


Ap(x) 


00 

P(a)e ^ ax da 

— oo 


then p = 0 for x < 0 implies that P(a) = P + (a), i.e. that 
P(a) is analytic in an upper half plane. But the relation 
between $(a), the Fourier transform of A<f>(x), and P + (a) is 
(for the case of uniform flow) 


$(a ) 


-j P + (a) 

r i , a M i 
P 0) [1 + — ] 


Since the denominator in the expression for $(a) has only a 
lower half plane zero (at a = - k/M), it is clear that $(a) is 
also analytic in an upper half plane (i.e. $(a) = $+(a)). 

This means that AtJ> will be zero for x < 0 and r = 1. We have 
verified that our solutions to the inlet and exhaust problems 
are compatible with other solutions for the case of uniform 
flow. 


III. DISCUSSION OF SAMPLE CALCULATIONS 


At the outset, let us clearly designate the convention 
we will be using in labeling the various modes. We shall 
always designate the lowest zero of a given lobe order m as 1 , 
i.e. n = 1 will be the lowest propagated mode of lobe number 
m. Then, a label such as, the (m, n) mode will mean an m 
lobed spinning wave and the n^* 1 radial order. (0, n) waves 
will be the symmetric waves, and the plane wave will be (0, 1). 
When we designate the phase of the reflected wave 0™ s , this 
will be taken to mean the reflected radial mode of an 

(m, n) incident wave. 

In the development of the theoretical model, a number of 
approximations were made, some physical and some mathematical. 
Each of these will clearly have their impact. Consider first 
the physical approximations: The most crucial physical 

approximation is the exclusion of any physical inhomogeneities 
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such as turbulence and neglect of jet spread and mixing. While 
this rendered the mathematics tractable, the use of the uniform 
plug flow model entailed the use of the particle displacement 
match boundary condition. Similarly, the thickness of the 
exhaust duct and the termination of the annular centerbody 
were excluded. The primary effect of these approximations is 
to render the phase calculations inaccurate, since the phase 
is particularly sensitive to such termination effects. The 
least effect of such approximations is on the amplitude of the 
radiation pattern. Thus it would appear that in phase synthesis 
of a radiation pattern composed of a set of flow duct modes is 
practically difficult if not rendered meaningless. On the 
other hand the restrictive effects of the mathematical approxi- 
mations introduced by the Carrier-Koiter procedure are probably 
much less (except near the axis as ® 0 for the radiation 

pattern, as discussed below) provided one carefully observes 
the rules devised by Koiter [9] and Carrier [10] . Indeed in 
our calculations this approximation is restricted only to the 
equivalent of the Wronskian normalized by the uniform flow 
Wronskian, so that for uniform flow and no flow, our theory is 
exact. In what follows, we shall try to illustrate these 
points by concrete comparisons, especially those involving 
experimental measurements. 


There are several interesting features displayed by the 
solutions such as equations (42a, b). The most obvious is the 
convective amplification due to the main jet, (1 + a M^/k^) 2 , 
and the weakening of that factor by the external flow, 

(1 + a M 2 /k 2 ) 2 . The second feature relates to the zone of 
relative silence defined by 


® < ® = cos-K 


/r^M| 


[M 0 + 


(1 - 


M 2 ) 


2 c 1 (1 + 1 ^) 


F 


M 2 


[M 0 + 


(1 - 


Mfy 


2 c 1 (1 + M 1 ) 


or for Mg = 0 


® = cos l { 


V C 1 

1 + M„ 


within the zone of relative silence the signal suffers attenua- 
tion due to the exponential factor_exp{Y-jy , n' arc cos(aQ-^/ic 2 ) ) , 
which results in a cusp at ® = ® s for the plane wave and 
nodes for the higher modes. The cusp has been previously 
observed by Gottlieb [16] . 


In the case of uniform flow (M-^ = M 2 ), the limit of the 
solution for pg as ® ■+ 0 is the same as the solution obtained 
for p^ (taking the case of the exhaust flow for concreteness) 
with the closure contour in the lower half a plane. Such is 
not the case for f Mg in the present solution. Indeed, the 
pressure distribution within the jet, as represented by, for 
example, (46) displays a zero pressure on the axis which_Ts 
not reproduced from the far field radiation pattern as ® + 

0. This appears to be the consequence of the Carrier-Koiter 
approximation. 
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Consider now the solutions for the symmetric modes, (0, 
n). For the case of no flow and plane wave incident, our solu- 
tion is exact as is illustrated by Figures 7, 8, 9. These re- 
sults may be compared with the results of Levine and Schwinger 
[2] among others. 

More useful comparisons, however, are offered by the 
experimental data of Plumblee et al. [6] and Mechel et al. 

[7] . Figures 10, 11 compare the flow and no flow data and 
includes the calculated effect of inclusion or exclusion of a 
hub. Plumblee 's data correspond to a monochromatic exhaust 
mode and a hub which is terminated in the duct exhaust plane, 
so that the calculated situation does not correspond exactly. 
The measurements have been made with great care. Not sur- 
prisingly, the no flow results show a reasonable agreement, 
however the cusped lobe of the calculation with flow is 
smoothed into an outward leaning cleaner lobe. Such a 
smoothing of the radiation pattern is to be expected in view 
of the scattering the radiated sound pattern by the intense 
turbulence of the spreading exhaust jet. This process is 
well illustrated by the calculations of Mani [17] , thus note 
in particular Figures 5, 6, 7 of that reference. Nevertheless, 
the influence of the exhaust jet refraction is well brought 
out by both the data and calculations. Indeed accounting for 
the diffusive smoothing of the pattern, the peak direction 
does seem to be approximated. On the other hand, the 
influence of the hub is seen to be negligible. 

Plumblee et al. [6] have also presented some measured 
radiation impedances. A sample comparison of these cases 
calculated on the basis of our theory are presented in Table 
I. The agreement, even in the case of no flow is somewhat 
less than satisfactory. This well illustrates the influence 
of the truncated hub in altering phase relations. Calculations 
without the hub included show no better agreement. What 
limited agreement there was with no flow, breaks down in the 
presence of flow. Thus the physical flow model appears to 
greatly exaggerate the acoustic mismatch between the duct and 
the surrounding, and the error, as will be seen presently, 
does not wholly originate of the Koiter approximation. 

Before considering the data of Mechel et al. [7], recall 
that the use of the uniform plug flow model resulted in the 
use of the Miles [12] - Ribner [13] particle displacement 
match boundary condition. In essence, this boundary condition 
results from the assumption of a clearly defined flow inter- 
facial streamline with the acoustic motion being transmitted 
from one region into another by the transverse perturbation 
of that clearly defined interfacial streamline. However, in 
the case of the ill-defined, highly turbulent flow of the jet 
boundary, the mean velocity changes "almost" smoothly across 
the interfacial region. In this case, noting the calculations 
of Mani [17] , there appears to be little reason to presume the 
picture based on the Miles-Ribner argument. Interest in this 
question was recently revived by Ingard and Singhal [18] who 
presented intriguing experimental evidence on the coupling of 
a sound source with a turbulent flow duct. Their evidence, 
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moreover seemed to favor the discarded notion of continuity 
of particle velocity rather than the displacement. 

It was hoped for that the present calculations be useful 
in this regard concerning the exhaust jet. The thought was 
that the duct reflection amplitude is a function of the boundary 
condition used. Then a comparison of the low frequency pre- 
diction using the two boundary conditions when compared with 
low frequency flow duct experiments of a similar nature might 
sort out this thorny question. Indeed the data of Mechel et 
al . [7] offer this opportunity, and it was a simple matter to 

modify our theory to include a particle velocity match at the 
jet interface. 

Such a test, as discussed above, is offered in Figure 12. 
Both the boundary conditions show a similar trend and predict 
reflection coefficients higher than the uniform flow case. 
Neither calculations reproduce the peak in measured reflection 
coefficients, but both do reproduce the reflection coefficients 
to within an acceptable error. However the marginal variation 
in the calculation does not allow a critical judgment of the 
boundary condition. As a closing remark for this aspect of 
the discussion, it should be noted that the no flow results of 
Mechel et al . did reproduce the calculated results to a high 
degree of accuracy. 

Finally, before proceeding to the discussion of asymmetric 
modes, the influence of the Mach number and flow mismatch are 
illustrated in Figure 13. 

Consider now the asymmetric modes. In light of the above 
discussion, we shall concentrate on results for the radiation 
patterns. In Figure 14, a comparison between our theory and 
the experimental results of the phased array model of Plumblee 
et al. [6] is offered. A disconcerting feature of these data 
is the non-zero radiation amplitude on the axis even without 
flow. This suggests spurious effects not accounted for in the 
physical model on which the experiments are based. Similarly, 
the theoretically predicted lobe in the shadow zone of the jet 
cannot be picked up over the flow noise of the acoustic pickup. 

A somewhat better agreement is apparent for the calculated 
(4, 1) mode in Figure 15. Plumblee 's data have not been 
plotted on this figure and reference must be made to Figure 20 
of that paper. However the refractive shifting of the pattern 
and the narrowing of the dominant lobe are both reproduced. 
Figure 15 also presents the effect of a weak external flow 
which may be interpreted also as a flight effect. The ultimate 
situation of the flight speed matching the jet speed being the 
uniform flow case. The effect of the weak external flow is to 
shift the pattern towards the axis and reduce the shadow zone 
lobe . 


More dramatic example of the jet Mach number effect is 
shown in Figure 16. One sees here the conflict between the 
axial convection at high Mach numbers, beaming or trapping 
the dominant component of the radiation in the shadow zone and 
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the sharp refraction towards the sideline. The flow enhanced 
frequency k 1 also increases the number of nodes. Compare this 
figure with the plane wave calculation in Figure 2 of Mani [5] . 

The physical model we have chosen is expected to better 
represent the exhaust conditions. However, the effect of flow 
inhomogeneities has been to weaken the correspondence. On the 
other hand in the inlet mode, the smoother sink flow is not 
represented by our model. Yet paradoxically, note the com- 
parison of our calculations with the Crigler and Copeland [8] 
data. It would appear that the 53 - 62 rotor/IGV interaction 
has resulted in a dominant (9, 2) mode. Compare too the 
phased array data of Pasko [19] with our calculations in 
Figure 18. In both the cases, we used a small external flow 
to weaken the degree of the flow mismatch. 


In the turbine exhaust noise problem, the dominant flow 
mismatch is the temperature. From simple kinematic con- 
siderations, it is possible to illustrate the influence of the 
hot jet is to sharply beam the sound into the sideline. Thus 
Candel [20] has shown that a plane wave incident on a flow/ 
temperature interface at an angle 0.. is sharply refracted to 
an angle 0 2 


cos 0£ 


(c 2 /c 1 )cos Q 1 
1 - cos 0^ 


with M 2 = 0. Hence consider the radiation patterns of Figures 
19, 20. Figure 20 is somewhat more meaningful since it cor- 
responds to a physical case of fixed axial velocity and fre- 
quency of the upstream sound source. 


IV. CONCLUSIONS 


In this study we have offered a method of calculation 
for spinning acoustic duct wave guide mode radiating through 
a flow mismatch. The convected wave equation was solved in 
part exactly, partly accounting for the asymmetry introduced 
by the flow mismatch into the mathematics by the Carrier-Koiter 
procedure. An attempt was made to compare the predictions of 
the current theory with realistic experimental models. It is 
concluded that the value of the calculations is primarily in 
explaining qualitative features of spinning mode radiation 
patterns. Real flow features not accounted for especially 
preclude the reproduction of phase information. Paradoxically 
there appears to be a better qualitative agreement for the 
inlet mode for which the physical model is not really meant 
than the exhaust mode which we tried to approximate. This 
appears to be the result of relatively cleaner inflow condi- 
tions. An attempt was also made to test the issue of the flow- 
acoustic boundary condition, but the result was not decisive. 
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Figure 1 Schematic of the model problem. 
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Figure 2 Branch cuts on the complex a -plane. 
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Figure 4 Asymptotic approximation of Q(x) function. 
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Figure 7 Plane mode reflection coefficients. 
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Figure 8 Plane mode reflection coefficients. 
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kj = 4.58 
FLOW = 150 FPS 



Figure 11 Comparison of a hubless plane wave radiation pattern 
with data of Plumblee [6], 



el et al , experimental data 
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Figure 12 Comparison of experimental and calculated reflection coefficients. 




Figure 14 Comparison of theoretical asymmetric (3, 1) radiation pattern 
with data of Plumblee [6]. 










Figure 17 Trial calculations for an inlet mode and comparison with Crigler 
and Copeland [8] fan data for 53 bladed fan and 62 bladed IGV (9 lobe 
spinning mode). 
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Figure 18 Annular inlet directivity patterns, comparison of theory 
with phase speaker array data of Pasko [19]. 
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Figure 20 Effect of hot exhaust jet with fixed physical exhaust jet speed 
and sound frequency. 
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